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Scatterings, particularly those involving resonances, and other elementary 

O 

■ processes do not happen instantaneously. In the context of semiclassical nu- 
O' 

I/"} ■ clear reaction simulations, we consider delays associated with an interaction 

ON i 

for incident quantum wave-packets. As a consequence, we express delays 
q ■ associated with elementary processes in terms of elements of the scattering 

matrix and phase shifts for elastic scattering. We show that, to within the 

■ second order in density, the simulation must account for delays in scattering 

■ consistently with the mean field in order to properly model thermodynamic 

properties such as pressure and free-energy density. The delays associated 

with nucleon-nucleon and pion-nucleon scattering in free space are analysed 

with their nontrivial energy dependence. Finally, an example of s-channel 

scattering of massless partons is studied, and scattering schemes in nuclear 

reaction simulations are investigated in the context of scattering delays. 
PACS numbers: 24.10.-i, 24.10. Cn, 25.75,+r, 21.65. +f 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



Semiclassical transport-simulations have provided the theoretical backbone for interpret- 
ing heavy-ion reactions at energies in excess of few tens of MeV per nucleon. Simulations 
follow either a cascade approach |T]j2|] or a solution to the Boltzmann equation including 
the mean field and the Pauli principle Within these approaches particles move clas- 

sically in-between collisions where momenta change abruptly and possibly new particles 
form. In practice, the Pauli principle appeared primarily important for reactions below 
150 MeV/nucleon. For beam energies of hundreds of MeV per nucleon, detailed compar- 
isons of simulations with data have been carried out in order to determine features 
of the momentum and density dependence of the nuclear mean field. Within many-body 
theory, efforts were made to calculate nuclear optical potentials starting from elementary 
NN interactions |T0|-[12|. In describing reactions up to the beam energies ~2 GeV/nucleon, 
pions and low-lying resonances have also been included in simulations. 

It has been noted that results of heavy-ion reaction simulations can depend sensitively 
on prescriptions of space-time details of elementary scattering processes |13|-|16|,P,|l7| . Such 
details of elementary processes may matter even more at higher beam energies than at 
those of the MSU or SIS accelerators, specifically, at those of AGS (14 GeV/nucleon), SPS 
(200 GeV/nucleon), and of the constructed and planned colliders, RHIC and LHC. This is 
because of high Lorentz dilation factors that can amplify space-time effects, and because of 
high particle densities early in the reactions. 

In this paper we investigate the time duration of elementary processes in free space, 
in the context of semiclassical simulations. Both the forward-going and scattered waves get 
delayed in elementary processes. While delay of the forward wave is consistent with effects 
of an optical potential on a particle in a medium, to the lowest order in particle density, 
both types of delay times, of the forward and scattered waves, affect the density of states in 
energy, to the lowest nontrivial second-order in particle density. As one important outcome 
of our work, we find that, to describe the density of states and associated thermodynamic 
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quantities such as pressure, heat capacity, etc., the delays for scattered particles must be 
properly accounted for in semiclassical simulations, in addition to the effects of the optical 
potential. If one relaxes the requirements on interactions in simulations, forsaking e.g. the 
proper behavior of single-particle energies, but demanding that free energy as a fundamental 
thermodynamic function, is properly reproduced to within the second order in density in 
equilibrium, then one can manipulate the mean field to entirely absorb effects of the delays 
for scattered particles. Alternatively, time delays for scattered particles may be manipulated 
to absorb the effects of optical potential. Even in a degenerate Fermi gas, as we show, the 
delays in scattering can affect thermodynamic quantities. We investigate several examples 
of particle-particle scattering amplitudes describing, respectively, the scattering through 
a sharp Breit-Wigner resonance, the irN interaction in the A channel, and the s-channel 
scattering of massless partons. Using free-space amplitudes, we calculate numerically the 
time delays for the NN system. We show that the delay for scattering has some dependence 
on the scattering angle. The delays become small at high energies. Finally, we investigate 
scattering schemes in simulations PIpiySJTB-H] in the context of time delays. The schemes 



include hard-sphere scattering and scattering at the distance of closest approach. We outline 
numerical strategies for determining the equations of state and transport coefficients to 
an arbitrary order in density given the scattering prescription. 

The content of subsequent sections is as follows. In Sec. || we consider a wave packet 
traversing a region with a scatterer and derive expressions for the delays of scattered and 
forward-going waves, in terms of scattering amplitudes (phase shifts) and their derivatives. 



In Sec. m we demonstrate that these delay times are consistent with an ergodic constraint 
stating that the extra time spent in the vicinity of a scatterer should be proportional to 
the change in the density of states brought about by that scatterer. We further show that 
the delay time for a forward-going wave corresponds to the classical motion of a particle in 
a mean field given by the self-energy for a plane wave in a low-density medium. We also 
discuss in that section and more in the Appendix, in the context of ergodic constraints, 
the modification of the scattering of two bodies brought about by the presence of a third 



body. Finally, we express the free-energy density and pressure in terms of different time 
delays. In Sec. [IV] we draw conclusions on the time delays in a degenerate medium by 
examining an expression for the change in the pressure due to two-particle interactions. 
Section [V] is devoted to the examples of resonance scattering. In Sec. [VT| we calculate spin- 
isospin averaged time delays for NN system and directly relate the delays to the second 
virial coefficient. Different scattering prescriptions and methods for the determination of 
equations of state and transport coefficients, to within arbitrary order in density for given 



prescriptions, are discussed in Sec. [VII 



II. TIME DELAYS FROM WAVE-PACKET DYNAMICS 

We consider a wave packet incident on a large spherical volume of radius R, at the center 
of which interaction takes place. To make the derivation simpler we assume that this volume 
is much larger than the wave packet. The wave function is 



Mt) = J dE g{E) VC(t), (2.1) 

where for large distances 

= E<^M) X (3 
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Here a and (3 denote the initial and final channel, respectively, and x is the internal wave 
function. Most often we shall consider elastic scattering of particles with no internal degrees 
of freedom, hence 

f(0, £) = -^E (2* + 1) Pt{<x*0) %(E), (2.4) 



where the amplitude % is related to scattering matrix and phase shift by 

Tt = i(Si — 1) = — 2sin<^expi<^, (2.5) 

Si = exp 2i5e. For the asymptotic form to be valid r must be much larger than L max /k and 
Tn must vanish for t > L max . 

The average time t vo i that a particle spends in the volume can be written as 

Tvol = ti — , (2.0) 

where current density in channel 7 is 

j 7 = 2^;(r,i)30 7 M), (2.7) 

and N is integrated incident flux through the surface. The time in (|2.6|) is represented as 
the average time for the particle to leave the volume minus the average time for the particle 
to enter the volume, as j • f is positive when a particle exits and negative when a particle 
enters. We decompose j into incoming and outgoing pieces and, as there is no interference 
between incoming and outgoing waves for a large volume, we write explicitly 

Tvol = T ou t Tin- (2-8) 

With the incoming wave having only a contribution from the unscattered wave packet, one 
finds, from (|2.3| ) and (|2.6|) , that r in = —R/v a . 

The outgoing wave has contributions from both the scattered and unscattered portions 
of the wave packet. Correspondingly, the outgoing current has three contributions: from 
the scattered wave, from the unscattered wave, and from the interference between the two. 
We first calculate the exit time associated with the scattered wave in the single-channel case, 
for particles leaving at the center of mass angle 6, using partial-wave expansion, 



r s {9) = R 



2 fdtj s (R,8,t)-rt 
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where Z is the same expression as the sum on the r.h.s. of the last equation, only without 
the time t in the integrand, 

Z= J2 dtdE dE> 9( E ) 9( E ') e i{k - k ' )R (2£ + 1) P^(cos 0) 
x(2£' + 1) P e {cos6) V(E') Ti(E) e -*C^>* Q + I) . 
By making the substitution, 



(2.10) 
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then carrying out integration by parts to make the derivatives with respect to E and E' 
act on g, expi(k — k')R, and T, and further integrating over the time coordinate and one of 
the energy variables, one obtains a simple expression for r s . For a well-defined energy this 
expression is 

t.(6) = - + At s (6) 
v 
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where 
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Z' = J2 ( 2 ^+ l)P*(cos0) (2/ + 1) P// (cos 9) T£(E)T t (E). 
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The more general result when many channels are open, in terms of the scattering ampli- 
tude /, is 
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r* dfg/3 _ df*p 



(2.14) 



v p 2i\f a tf { Jap dE J ^ dE , 

This is of use when determining time delays for a NN system in Sec. [VT|. In that case the 
indices refer to spin components. 

For the case of elastic scattering in only one partial wave, the exit time from ( |2.12p may 
be written as 



Thus, the extra time delay due to the scattering is dbjdE which is not the naive guess 
that one would make from considering an incoming partial wave reflecting off a potential. 
In that case the answer should be 2 d5/dE since the potential modifies the scattered wave by 
a factor e 2lS . This discrepancy is associated with the interference between the forward-going 
wave and the scattered wave. 

Next, we consider the current and the average exit time associated with the interference 
of the forward wave 0° with itself and with the scattered wave. The angular integration 
must extend up to 9 ~ L max /kR which limits the forward wave, and it may, in particular, 
extend over the whole angular range. Analogous procedure to that before yields 

D 1 Lmax A 

x / = --^ 7 g(2. + i)-(r«( B ) + r; (B )) (2.16) 

where 

Lmax 

Z" = £ (2£ + 1) (1 - i{%{E) - %*{E))) (2.17) 

In terms of phase shifts this gives 

= R Ttl (2* +1)2 00s 2ft § 
Tf v YJi max (2i + 1) (1 - 4 sin 2 S e ) ' 1 ' ' 

The more general result in terms of the scattering amplitude / is 

t? = — + —r ; ^ —-^lk a Ref aa (0)}. (2.19) 

By combining the results for scattered and forward waves, we obtain the average exit 
time in the single-channel case equal to 

N s \ r , n .dN a /dn 



Tout = Tf(l- j^J + J d£lT s (tt) 



N 

R S, W (2l+l)2ff 
v J2e max (2£+l) ' 

The total number of scattered particles is iV s . The delay time ( j2.2C|) is the same that one 

would have guessed assuming that different partial waves acted independently and each was 



delayed by 2d5e/dE. The complexity of Eqs. ( |2.18|) and ( |2.12|) stems from interference of 
partial waves. In the next section we discuss the consistency of these results with equilibrium 
expectations. We show that the forward delay time agrees with a delay in the motion through 
a mean field. 



III. THE ERGODIC CONSTRAINT AND THE MEAN FIELD 

A. Ergodicity 

All thermodynamic variables such as the pressure can be found if one knows the density 
of states within a system. Assuming that particles do not interact, and that statistical effects 
can be ignored, one can obtain thermodynamic quantities to lowest order in particle density. 
Assuming that particles interact only two at a time allows one to calculate these quantities 
to the next order in the density. To account for the interaction of the two particles at a time, 
one needs to find the correction to the density of states of relative motion. This correction 
to the density of states p{E) = dn/dE is given by phase shifts for two-particle scattering 



p(E) = Po (E) + Ap(E) = + 1 £ (21 + 1) §, (3.1) 

(27i) d v Ti ^ dE 

where k is the reduced relative momentum, V is the volume, and v is the velocity of relative 
motion. This is a complete quantum-mechanical answer in the single-channel case. If pro- 
cesses such as a + b — > c + d can occur, a similar answer can be obtained by diagonalizing 



the S-matrix and finding the eigenphases [20 



In a simulation, the equations of motion should result in a change in the probability for 
two particles of a specific energy to be in the vicinity of each other, proportional to the change 
in the density of states for these particles. If the equations of motion accomplish this, the 
classical simulation will yield correct thermodynamic quantities. We consider a subvolume V 
in relative-coordinate space and a narrow relative-energy range \E — Eq\ < 5E. Given that 
particles spend a portion of time tq within V in the absence of interactions, when sampling 
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is carried out over a long time T, then, in the presence of interactions, the time within V 
should change by At such that 

Ar = Ap = I S< (2l + l)g 

TO PO PO 

This can be thought of as an ergodic constraint, since a system continuously sampled over 
a long time should be populated according to its contribution to the density of states. In this 
section we explore the case of elastic scattering and wish to see if the results of the preceeding 
section are consistent with the ergodic constraint and with the picture of particles moving 
through a mean field. 

In the context of a classical simulation, the additional time Ar will come from three 
causes. First, particle kinetic energy within the volume V will be different as compared to 
free space, by the negative of the mean potential within the volume, —u(k). The change 
in the energy gives a change in the momentum and in the velocity and, in consequence, 
in the time spent within subvolume by an amount denoted &.r{ las . In addition, the en- 
ergy dependence of the mean potential generally changes the velocity by Av = du/dk. 
The corresponding change in the time spent within the volume is denoted as Ar| /as . Both 
these contributions to the time follow from classical equations of motion within a potential. 
A third contribution At coU stems from collisions. Because of collisions the particles will 
emerge on the average earlier or later (particularly in the case of resonance scattering) from 
the volume than in absence of collisions. Summing the three contributions should yield 

Ar clas + Ar clas + ^coll = 1 E , ( 2 f + 1) § 
TO PO 

The time Arf as can be easily obtained from the ergodic theorem for a classical potential. 
Otherwise, one can resort to geometric considerations caring, however, for the fact that, 
to the lowest order in u, the velocity changes its direction as well as the magnitude. As the 
relative change in time spent within the subvolume should be proportional to the relative 
change in the density of states brought about by the change in kinetic energy, we get 

Arf- = Apo = M£)ff 

TO PO PO 
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The contribution Ar| /as is obtained by noting that the time it takes to traverse a given path 
through the subvolume is proportional to the inverse of velocity, 

A das Ay I du du > v 

= = — 7i: = -7^- ( 3 - 5 ) 

To v v dk dhj 

The contribution due to collisions, At co11 , will be determined by the collision rate and the 
change Ar s in the average time spent within the collision range compared to free flight time, 

Ar coll = r yAr s . (3.6) 

Combining the expressions for different contributions to the change in time on account of 
interactions, we find that, according to the ergodic constraint, we should have 

n ? (2£ +1) dE = -dE U - po dE +Po V (3 - 7) 
In the following we shall express u in terms of phase shifts and then determine At s 
from ( |3.7| ) and examine whether the result is consistent with what was obtained in preceeding 
section. The potential u in (|3.7|) is identified with the correction to the kinetic energy at 
which the single-particle Green's function has a pole, 

1 



Ei -e(ki) -w + 27/2 

' (3-8) 



- 1 



where we use the fact that self-energy can be expressed in terms of the forward element of 
the T-matrix and where k continues to be the relative particle momentum. Assuming that 
angular momentum is conserved, we obtain 

(k|T|k) = ^ E ( k K m ) ( £ > m \ r \ £ > m ) ( £ > m l k ) 

= E(^ + 1 )^T- (3-9) 

The imaginary part of the matrix times —2 can be shown to be equal to av by using the 
standard expression for the cross section involving sin 2 6. Then 7 is indeed the geomet- 
ric scattering rate one expects. Turning now to the real part of T-matrix, the potential 
associated with the presence of other particle within subvolume becomes 
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SAM +!)»■■» 2ft (3 . 10) 
27rp 

The sum of the times 

AT das and Ar das in Q invo i ves t h e derivative d (u(E)p (E)) /dE. 

With ( EHUD , we obtain 

A das I \ T clas 1 J 1 

1 + 2 = " 1 ^ (Po tt(^)) = — E (2* + 1) cos 25, §. (3.11) 



Upon inserting the above result into Eq. fl3.7|) we find 

_ £,(21+1) sin^g 

E£ (2£+l)sin 2 5 £ • ld - i2J 

If only one partial wave is scattered, then the change in time reduces to 

d5t 

s = ii- < 313) 



This result is consistent with Eq. ( |2.15| ) of the last section. If many waves are scattered, 



then there is generally interference between various scattered partial waves giving rise to 
a variation of the change in time with scattering angle, cf. Eq. ( |2.12| ). Equation ( |3.12| ) 



represents, in such a case, an average of the change of time over the scattering angles, 
weighted with the flux or cross section. 

One can further see that the time delay of the forward wave, as derived in the last section, 
may be identified with Arf as +Ar 2 das . The wave packet in the last section was of a transverse 
size B = L max /k. We further assume that the length of the packet is W. The time for such 
a packet to move by the potential in the absence of interactions is tq = W/v. Using Eq. ([2.1 6|) 
in the limit of large L max we obtain 

Ar ; gg (2* +1)2 cos 25, § 

?5M2f+l) ' 1 j 

For wide packets we can approximate J2e (2€ + 1) ~ 2 J idl = k 2 B 2 and, with the volume 

of a packet being equal to V — nB 2 W, we can rewrite ( |3.14|) as 



Ar f _ Ee (21 + 1) cos 26 1 § _ E £ (21 + 1) cos 25 e § 

where p is the density of states inside V. One can see that this agrees with (|3.11|) . 
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Given a uniform many-body system, with density n and temperature T treated as inde- 
pendent variables, all intensive thermodynamic quantitates can be obtained from the free 
energy per unit volume f . In terms of the correction to the density of states of relative 
motion in energy in ( |3.1| ), the free energy f, to the second order in the density, can be 
written as 

f(n,T)=f (n,T)+AP(n,T), (3.16) 
where fo is free energy for a noninteracting system and AP is the correction to pressure 

P = * + M = »T-T* (^f\j iEe -^ {2e+l) ^ (3 ,T, 

and where the effects of statistics are ignored. It can verified by a direct calculation that 
the correction to pressure may be further expressed as 
1 f dpi dp2 



Ap = 2 / Wf /(Pl) Re ((Pl " P2,/21 7 l(pi " P2>/2> 

- T \lwfWf !M! ^l itti ^ T - (a) - (318) 



Here / = E ^ T is phase-space occupancy. The first term on the r.h.s. of ( J3.18 ) accounts 



for the forward time delay or mean field. The second term accounts for delay in scattering 
and, depending on the sign of Ar s , allows further for an interpretation in terms of the 
reduction in the number of degrees of freedom or in terms of excluded volume. 



B. Limitation of the Considerations 

Two important conclusions can be reached regarding the method of the previous section. 
First, the sum of time delays for scattered and forward-going waves is consistent with ergodic 
constraints. Secondly, the time delay of the forward-going wave is consistent with motion 
of a particle through a mean field. However, this consistency was derived assuming that 
the size of the region used to determine the densities was large. The limitations of this 
approximation are discussed below. 
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A problem with incorporating a mean field which depends on the number of particles 
within a given volume element, is that such field can result in classically bound states for at- 
tractive potentials, or inpenetrable potential barriers for repulsive potentials. Considerations 
in this section required that the potential u(E) was much smaller than any characteristic en- 
ergy. This is necessary for Ap <C p . For a sufficiently large volume V this is not a problem. 
However, for a finite V there exists phase space that is either bound (attractive potentials) 
or unavailable (repulsive potentials), of a magnitude 

A = j^f (2mu(k = 0)) 3/2 . (3.19) 

On using sin^o ~ —ka, where a is scattering length, for low k in the expression for optical 
potential (|3.10|) , one finds that 

4 a 3 / 2 

Thus one finds that for a strict validity of our considerations the volume V containing a single 



scatterer must be chosen large compared to the scattering length. On examining Eq. Q2.16 ) 



from the previous section, one finds that, in the first place, such a condition must be satisfied 
to allow for a proper definition of Tf. Scattering lengths are of the order of the interaction 
range, unless a resonance or bound state exist close to the threshold. Unfortunately, this 
is precisely the case for nucleon-nucleon scattering near threshold where scattering lengths 
approach 20 fm. 

C. Scattering in the Presence of Third Bodies 

When considering the interaction of two particles in Sees. H and |III A| , we neglected 
the probability that a third body could interact with any of the two particles within the 
subvolume used to define the mean field. In Eq. Q3.3p the contributions to the extra time 
spent within the subvolume, Arf Zas and Ar|' as , arose because the particle trajectories were 
modified by the mean field. The first contribution arose because the kinetic energy and 
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direction were altered on entering the mean field. The second contribution arose due to 
dependence of the mean field on momentum that led to the change of velocity by du(E)/dE. 
When many scatterings occur, the time spent in a given relative state does not just depend 
on trajectories and velocity, but further on the manner in which such state is populated and 
depopulated. If population of different states is to be consistent with first two of the time 
delays in Eq. (|3.3|) due to interactions, the scattering prescriptions need to be modified in 
the presence of third bodies. 

Transition rates per unit time from an initial state, are given though, generally, by a tran- 
sition matrix element to a final state squared multiplied by the final-state density in energy. 
For example, in the case of particles 1 and 3 scattering in the presence of particle 2, the final- 
state density would be obtained from a product of single-particle densities for 1 and 3. Latter 
densities, up to a factor, are identical with imaginary parts of single-particle Green's func- 
tions ( |3.8|) . When shifting the single-particle energies by mean fields generally dependent 
on relative momenta of 1 and 2, and 3 and 2, one precisely accounts for the change in the 
density of relative states, in feeding of the states. Besides the change of relative momenta 
due to mean field, that corresponds to /S.rf ass , the levels of 1 and 2 or 3 and 2 are pulled 
apart or pushed together when mean field is energy- dependent, changing the density, what 
corresponds to Ar| Zass . Then the rates calculated with single-particle energies in the fi- 
nal density of states, would yield population of states consistent with these first two times 
in Eq. (|3.3|) , if the particles could leave these states during the whole time such as obtained 
with an inclusion of the forward time delay or simply adoption of mean field. 

Ergodicity when many scatterings occur is discussed in more detail in the Appendix, 
together with the single-particle spectral densities. We show there, that, in fact, population 
of different states could be made consistent with an ergodic constraint involving all three 
delay times in (|3T3|), i.e. also the scattering delay, if one included in the Green's function fl3.8| ) 
the imaginary part of the T-matrix and allowed for a separate dependence of the matrix 
on energy and relative momentum. Modification of the scattering or transition rates in 
terms of such complete spectral functions, nonetheless, might not be practical except for 
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fully equilibrium situations. Also, it needs to be stated that a modification of the final- 
state density for two particles such as 1 and 3 on account of scattering with third particle, 
independent of a modification of matrix element for 1 and 3, might not be proper. With 
regard to scattering, correlations could be important, persisting throughout the interaction 
process. 

It will be seen in the examples of Sees. |V| and [VT| that the contribution to the density 
of states from the mean field alone can be very significant. If scattering is then done in the 
presence of strong mean fields, it is important that matrix elements get modified too. In case 
of momentum- independent mean fields, this just amounts to the calculation of scattering 
rates in terms of kinetic particle energies only. Situation can be more cumbersome in case 
of momentum- dependent fields as one generally loses the convenience of working in the rest 
frame of scattering particles. A particularly strong momentum-dependence of the mean field 
results for pions in nuclear matter due the derivative coupling to nucleons p2| , p3| . 

Rates and scattering prescriptions may be complicated when the time delay is negative. 
In such a case the prescription may effectively exclude a relative volume, e.g. an energy- 
dependent hard core. One may not want to fill the excluded phase-space volume in the 
scattering with a third body and with this populate relative states following ergodic theorem 



including the scattering delay time. Scattering prescriptions are investigated in Sec. [VTJ 



including one with a hard core prescription and one where effective time delays are generated 
by correlating the outgoing scattering angle with the impact parameter. 



D. Manipulating Time Delays 

If one solely aims at satisfying the ergodic constraint, and thus properly reproducing 
thermodynamic properties of a system, while ignoring the physical differences behind the 
mean field or forward time delay and the delays for scattered waves, then one can absorb 
all interaction effects, in the lowest nontrivial order in density, exclusively into the mean 
field or, alternatively, into the delays for scattered particles. We shall illustrate our points 
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by considering the correction to the pressure in fl3.17| ) and ( |3.18|) . 

Thus, from Q3.2D , ( |3.10| ) and ( |3.11| ), it follows that the density of two-particle states in 
energy would be properly reproduced when ignoring the delays for scattered waves and, 
in place of the mean field u in the center expression in ( |3.10| ) , using a field v! given by ( |3.10| ) 



with sin 25i replaced by 2(6t(E) — <5^(0)). As far as the correction to the pressure is concerned, 
this corresponds to rewriting the expression ( |3.18| ) as 

AP = \ J rfpi dp 2 /(pO f(p 2 ) (( Pl - p 2 )/2| Re T |( Pl - p 2 )/2) , (3.21) 

where 

<p|ReT» = ^( 2 £+ 1)(5,(0) - 5 e (E)). (3.22) 

The replacement of the sine by its argument in the field is actually a good approximation 
when the phase shifts are low compared to tt/4. The general conclusion then is that the 
delays for scattered waves are relatively unimportant for the thermodynamic properties when 
the phase shifts are low, no matter what are the values of these delays. 

If the mean field were to be ignored, then, in order to properly reproduce the density of 
two-particle states, from (|3.2p and (|3.6p, the mean delay time for scattering should be taken 



equal to 



s ~ 2£,(2£+l)sin 2 V 



and for scattering in only one partial wave 



A . 1 d5e 

A < = V~^T^ 3.24 
2 sin bp ah 



in place of ( |3.12| ). The correction to the pressure ( p.l^ ) is then rewritten as 

AP = -T l - J d Pl dp 2 f{ Pl ) f(p 2 ) a v Ar' s . (3.25) 

Some advantage of an approach with delays put into scattering is that particles move with 
uniform velocities. One problem, though, is that the delay times At' s diverge at thresholds, 
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where phase shifts go to zero linearly in the momentum with the coefficient of proportionality 
being equal to the negative of scattering length. In Sec. |V| we shall discuss resonance 
scattering including, of interest for heavy-ion collisions, pion-nucleon scattering, and in that 
context discuss more the different possible prescriptions for time delays. Now we turn to to 
the role of statistics. 



IV. TIME DELAYS AND STATISTICS 

For any finite energy the density of states in energy is increases by symmetrization and 
decreased by antisymmetrization. Corrections to thermodynamic quantities, such as the 
free energy density or pressure, arise within the second order in density even in the absence 
of interactions. Symmetrization also affects scattering processes both in that the outgoing 
states may be Pauli-blocked or Bose-enhanced and in that the scattering amplitude may be 



internally modified J24]. In the following, we consider pressure in an equilibrated many-body 
system in terms of the T-matrix. In a nonequilibrium system, for low scattering rates, the T- 
matrix approximation |25| in the single-particle equations of motion leads to the Boltzmann 
equation with rates corresponding to two-particle collisions, enhanced or reduced on account 
of statistics of final states and with medium-modified amplitudes. We assess implications 
for time delays and scattering processes in simulations, following from ergodicity. 

The pressure in a many-body system at a given temperature T and chemical potential fi 



may be generally represented as 



P(ji,,T) = P (ji,T)+AP(ji l ,T), (4.1) 

where 

dP dp Jrn 1 



» 1 f dP dr> 



x 



(27r) 3 (27r) 3 1 e {E t -2^/T _ 1 

jT 1 d\ (-i) Im (p|V G A (P, Et) (|p> ± I - p» . (4.2) 



Here V is the interaction, p and P are relative and total momenta, respectively, E t is total en- 
ergy, and G\ is the two-particle Green's function within a system with the interaction scaled 
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down by a factor of A. The negative of the imaginary part of this function, divided by it, 
plays a role of the two-particle spectral function. The statistical factor in (|4.2|) is bosonic, 
as is appropriate for a state of two particles with the same statistics. The apparent change 
in sign of the correction to the pressure above in the classical limit (fi/T — > — oo), compared 
to ( p.l7| ), is associated with the fact that there the pressure is expressed as a function of 
density while in (|4.2| ) it is expressed as a function of /j,. 

In the T-matrix approximation |25| the two-particle Green's function satisfies 



G = G + G VG, (4.3) 

and explicitly 

(k|G(P,^|ki) = (k|G (P,^)|k 1 )+/^^(k|G (P,^)|k 2 ) 

x(k 2 |V|k 3 )(k 3 |G(P,£; t )|k 1 ). (4.4) 

The noninteracting Green's function in the above is equal to 

(p|G Q (P, E t )\p') = (27r) 3 5(p - p') p2 !! (P,P) 2/ —■ • (4-5) 

h t — P z /4m — jr/m + le 

The factor of N stems from an equal-time commutator of the operators for two particles, 

iV(k, P) = (1 ± /(P/2 + k))(l ± /(P/2 - k)) - /(P/2 + k) /(P/2 - k) 

= l±/(P/2 + k)±/(P/2-k), (4.6) 

and the upper signs refer to bosons and lower to fermions. The T-matrix and Green's 
function are related with 

T = V + VGV, (4.7) 

and thus the T-matrix satisfies 

T = V + VG Q T. (4.8) 

The factor of N in Go in the above equation, on considering the case of a resonance, can be 
related to the fact that, for a bosonic two-particle state, the width is equal to the difference 
between the decay and the formation rates ||25|| . 



On integrating in ( f4.2| ), we obtain, from ( |4.3|) , 
/ dAIm(p|VG A (|p)±|-p)) =Im(p|^-(VG ) n (ip)±|-p)) 



n=l 

-Im(p| log(l-VG )(|p)±|-p)) 



(p|atan(V(l - VReGo) _1 ImG ) (|p) ± | - p)) 
(p|atan(7eimG )(|p)±|-p)), (4.9) 



where 1Z satisfies 



TZ = V + VReG TZ. (4.10) 

On introducing then the matrix 1Z that is hermitian in the spherical angle and may be 
diagonalized, 

n( P , n, n', p) = sga (E t - 2^) \n( p , p, n) n( p , p, q')\ 1/2 ( p , n\n(p, E t )\ P , n') 

= ]T r n { P , p) y n { P , p n) yjj>, p n'), (4.11) 

n 

where E t = P 2 /4m + P 2 /m, and {y n } form an orthonormal set in spherical angle (with 
a definite symmetry under inversion), we can rewrite (f4.2|) as 



APlll < ^'llwjif e<E+PV4 m -2,)/T _ i E P), (4.12) 

where <5„ = atan (— P 2 r n /8n 2 v). 

In terms of the phase shifts 5 n introduced above, the on-shell T-matrix can be ex- 
pressed as 

(p, fi|T(P, 3t)|p, JT) = sgn(P t - 2p,) |iV(p, P, fi) iV(p, P, fi')P 1/2 

x J2 sin 5„ e l(5 " y n ( P , P, fi) ^(p, P (4.13) 



On carrying partial integrations in ( 4.12 ), the correction to the pressure may be decomposed 



into the mean-field and scattering contributions that take the form, in terms of the T-matrix, 
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AP(/.,T) = -i/^ J ^/(pi)/(p 2 )Re(p|r(P, e ( Pl ) + e (p 2 ))(|p}±i-p}) 



1 r dP P 2 r dp f/,. ,,. > 



J 7^)3 J d^(f(pi)f(P2)^vAr s (l±f(p[))(l±f(p' 2 )) 



-/(Pi) /(pa) ^ «At, /(pi) /(p£)V (4.14) 



2 7 (2tt) 3 6m 7 (2? 
(icx 

where / = ( e ( e (p)-^/ T q= l)" 1 , da/dfi = (m/47r) 2 | (p|T(P, E t ) (|p') ± | - p')) | 2 , integration 
is over the spherical angle of 2tt, and the time delay for scattering is 

(^"a(7^) Im(loE<p|r(P - £ ' )(|p ' >±| - p ' >)) - (4 - 15) 

Following ( 4.1|) and (|4.14{) , in order to produce proper changes in the pressure in simula- 
tions, associated with two-particle interactions, it is necessary to include the effects of mean 
field on particle motion, and to delay the collision processes, according to the expressions in 
terms of T-matrix, with an internal symmetrization and symmetrization with other particles 
in the medium accounted for in the final and intermediate states (Eqs. ( |4.8| ), ( f4.5| ), ( |4.14[ ), 
and ( |4.15| )). Besides delaying direct collisions, where two particles alter their momenta and 
spins, by Ar s , it is necessary to delay exchange collision processes, where pairs of particles 
meet and interchange pairwise quantum numbers, by — Ar s . (The latter holds for particles 
of one statistics. In the case of particles of opposite statistics, the exchange collisions should 
be delayed by Ar s .) The exchange collisions may be thought of as processes |26| where two 
particles collide and, while in an intermediate 2p-2h state, encounter two more particles, 
ending up in the interchange of the quantum numbers. The cross section, from ( |4.14|) , is the 
same as for the direct collisions. Notably, there is no room for the exchange collisions in the 
Boltzmann equation that ignores the duration of interactions, since these processes leave 
the occupations of single-particle states unaltered. These processes should, nonetheless, ap- 
pear in a possible quantum Enskog equation, since their duration affects thermodynamic 
quantities. 



V. RESONANCE SCATTERING 
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A. Sharp Breit-Wigner Resonance 



An obvious example, that can serve to illustrate the time delays and different prescrip- 
tions, is that involving a sharp resonance. If the width of the resonance T is small compared 
to its energy Er, then the phase shift in the vicinity of Er is given by 

172 

t ^ 5 = ~E~E~R (5 ' 1} 

This yields the following time-delays for the scattered and forward-going waves, respectively, 

A d5 _ T 

s= dE = 2[(E-E R ) 2 + T 2 /4] ' ( j 

and 

f x ~o~ = ~ s = rp-4 a + rv4] ' (5 ' 3) 

where Ary is the time delay of the forward going wave packet with total cross-sectional 
area irB 2 = (vr/fc 2 ) J2e{^ + !)• Figure |l] illustrates the different times. It is seen that the 
delay for the forward-going wave turns negative in the vicinity of the resonance. For light in 
a dielectric medium |27] this corresponds to the increase in group velocity for packets with 
resonant frequencies. 

The energy-averaged delay for the forward-going wave Ary is zero. (This may be expected 
whenever phase-shift variation is limited to a narrow range in energy.) Weighted with the 
cross section, the average delay time for scattered waves over energy is equal to the inverse 
width coinciding with naive expectations. The actual delay time for scattered waves ( |5.2|) 



is twice as high at the resonance, see also |28[ , and it drops rapidly with energy when going 
away from the vicinity of that resonance. 

Time delay equal to the inverse width is further obtained when putting all delays into 
scattering, from ( [5.2D and ( |5.3|) , 

Ari = I . (5.4) 
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It is now a common prescription in relativistic simulations to use an energy-independent 
delay given by the inverse width. Provided that particles are not also propagated through 
a mean field and resonances are indeed of a Breit-Wigner form, such a prescription would 
yield correct thermodynamic functions up to the second order in density 



B. 7riVA System 

Different resonances in a nN system may be identified as independent particles. For the 
lowest P33 A-resonance, the width is comparable to the energy of the resonance above the 
ttN threshold and, as such, this width exhibits a significant energy variation. Interacting 
system of pions, nucleons, and deltas, is of considerable interest for heavy-ion collisions in 
the beam-energy range from few hundred MeV/nucleon to few GeV/nucleon. 

In this section we address three issues. First, by studying 7riV phase shifts we calculate 
time delays and indicate differences with the Breit-Wigner case above. Second, we point out 
that a quantum decomposition of the effective increase of the density of states represented 
by the time delay would include contributions from both the A and ttN components. Finally, 
we point out a practical difficulty in separating the A and ttN components. 

The time delays can be calculated given the experimentally determined ttN phase shifts 
|29| , |30|1 . The delay of the forward wave (divided by the fraction of particles that scatter as 
in ( |5.3|) ) and the delay of the scattered wave are shown for a 7f + p system in Fig. 0. In many 
prescriptions, time delays are only incorporated into scattering events, which would mean 
that the combination of the two contributions At' s is appropriate. One should note that this 
combination differs significantly from the Breit-Wigner result and is extremely divergent 
at threshold due to a rapidly declining cross section, cf. Sec. [Ill D| . If the mean field is 



consistently incorporated into a simulation, the appropriate delay is At s which behaves as 
the derivative of the phase shift with respect to energy which is approximately the scattering 
length divided by the velocity in the vicinity of the threshold. This threshold divergence is 
much weaker than that of At'. 
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The imaginary part of a A Green's function in a thermal system gives the number of states 



for the resonance per unit energy and volume, Ap" = —Img^/ir [ pT| , |25| . The requirement 
then of a consistency with the number of colliding ttN pairs gives the time At" during which 
the ttN pair should turn into the resonant particle, 

Ap"{E) a{E)vAr' s ' 



p (E) V 



(5.5) 



with a result that is the inverse of resonance width 32,33 



^=rm- (5 - 6) 

This time differs from At' s , such as in Eq. fl3.24|) , that would, in particular, involve the energy 
derivative of T; likewise Ap" differs from Ap. Interestingly, within any single spin-isospin 
channel, either correction to the density of states, Ap or Ap", integrates over energy to unity, 
i.e. one net state is gained. For sharp resonances, with couplings and widths independent of 
energy, there is no difference between At' s and At", as is apparent from Sec. |V A| , and there 
is no difference between Ap and Ap". 

On studying the density of states for pions or nucleons as in the Appendix, one finds 
a change, per nN pair, that is equal to Ap(E) — Ap"(E), where Ap = (l/ir)(d5/dE) in 
any single channel, i.e. one that precisely compensates the discrepancy above. Same type of 
discrepancy and compensation is found when applying the considerations to a system of pions 
and rho mesons, where 7r + 7r <-> p reactions take place. If one were to ask about, in the last 
system, how many p mesons decay into dilepton pairs, the answer would involve the density of 
rho states within the system, rather than the overall change in the density of states associated 
with the resonance formation, important for thermodynamic considerations. For certain 
questions one has to keep in mind that the time delays derived before correspond to the 
change in the overall density of states and not necessarily to the existence of the resonant 
particles. 

The time for the conversion into a resonant particle in scattering Q5.6| ) may diverge 
strongly when threshold is approached, which parallels the situation when the overall delay 
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time associated with the interaction is forced onto the scattering. (The conversion time for 
a spherical wave, that should be identified as Ar^ = 2nAp" oc T(E), tends to zero at the 
threshold.) In practice, manipulations of the conversion time, dividing this time between 
the scattered and forward waves, may pose more difficulty than the manipulations of the 
overall delay time, as negative conversion times cannot be simulated. 



C. s-Channel Scattering of Massless Partons 

An example, where time delays associated with interaction are relevant, is the collision 
of partons in ultrarelativistic nuclear reactions. Partons are copiously produced early on 
in reactions and the goal of simulating partonic cascades [j34] -|36|l is to determine the equi- 
libration time scale and initial equilibrated energy density. Most partons at midrapidity 
are produced far off shell and decay via bremsstrahlung. Thus, quantum considerations are 
necessary to establish the duration of processes. 

While we shall not tackle the general problem here, we shall try to gain insight by 
considering the simplified example of elastic scattering of two massless partons in the s- 
channel, assuming that the intermediate particle is massless. This turns out to be quite 
similar to the case of ttN scattering close to the threshold. Since there are no energy scales, 
to lowest order in perturbation theory the phase shift must depend only on the coupling 
constant a, 

tan S = —ca. (5.7) 



Then, the overall time delay for a spherical wave from (|5.7|) is zero. This peculiar result will 
arise from any theory with no energy scale since phase shifts are dimensionless. 

However, QCD acquires a scale A through renormalization that gives the coupling con- 
stant an energy dependence |3T|j: 



47T 

tan5 = Ca(E) = C /3 log(£7A 2 )' (5 ' 8) 



With this, the time delay for spherical waves becomes 
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dS 3 . 2 . 3 c 2 a(EY 

= Sin = — x r 

dE arE axE 1 + c 2 a(Ey 



Ar sph = 2 — = — — sin 5 = - ; 9 /r7n9 • (5.9) 



If all the time delay is put into the scattered wave, as in Sec. [Ill D| , then the correct time 
delay for the scattered wave acquires a particularly simple form 

A ^4^E- < 5 - 10 > 

Note that the time delay does not involve the coupling constant and behaves as inverse 
energy. If one chose the time delay equal to the inverse width for the intermediate state, 
then would certainly have obtained a time proportional to l/(aE), as the width would be 
proportional to a. 

Since an s-channel scattering involves an intermediate state very far from being on- 
shell, the questions involving the time delay for such a process may not be so crucial since 
these processes are rather rapid. Of a greater concern is the formation of partons through 
bremsstrahlung involving intermediate states which are nearly on-shell. Since such processes 
create the majority of soft particles in an ultrarelativistic pp collision, the issue of when and 
where such particles appear can greatly affect estimates of the initial thermalized energy 
density. Unfortunately, such two to three or more particle processes are outside the scope of 
this analysis, but similar problems have been addressed in the context of decaying hadronic 
resonances |38[ . 



VI. NUCLEON-NUCLEON INTERACTION 

The nucleon-nucleon system is one for which the quantum-mechanical scattering- 
amplitudes have been most carefully measured in physics. It represents the most relevant 
case of scattering for heavy-ion physics, where semiclassical simulations utilizing single- 
nucleon degrees of freedom are commonly used to model heavy-ion reactions. At moderate 
densities and high temperatures, when cluster formation is unlikely, it may be reasonable 
to assume that nucleons interact two at a time. Using phase shifts inferred from scatter- 
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ing data, one can calculate quantum mechanical scattering amplitudes and determine the 
appropriate time delays following the prescriptions outlined in the previous sections. 

First we present the calculation of the time delay time for the scattered wave at a given 
angle, At s (8), when averaging over the spin and isospin directions. Due to antisymmetriza- 
tion and the conservation of angular momentum, isospin and parity, the magnitude of the 
net nucleon spin is conserved in NN interactions. The S = amplitude for a given isospin 
is given by 

f=o {9) = * £ (2£ + i) P/(co8 0) % , (6.1) 

where only even or odd values of I are included. That is compensated by the factor before 
the sum in (|6.1| ) being twice as large as in fl2.4|) . Given the conservation laws, the S = 1 
amplitude is of the general form, 

f^\0, 0) = ~E {101 Aw) (* iji - A 1 A*'|J» <t>) H*j , (6-2) 

e 

where Tut'j = i(Su'j — and for the coupled waves diagonal matrix elements are Su = 
cos 2ej exp 2i5je, and off-diagonal elements Swj = i sin 2ej exp i(5je + 5j£>) , where e is mixing 
parameter. In carrying calculations up to a laboratory energy of Ei^ = 400 MeV, we use 
all partial waves with both i and j less than 5. The mixing of waves is not very strong in 
this region. The phase shifts are generated using a potential model that had been carefully 
fitted to describe the NN data p| . 

Using Eq. (|2.14j) , one obtains a delay time At s that depends on the scattering angle. 
This time is shown by different lines in Fig. ^ for several values of E^. Due to the averaging 
over initial spin and isospin and the amplitude antisymmetrization, the time is symmetric 
with respect to 90°. The time, averaged over angles, is additionally shown as a function 
of the laboratory energy by a solid line in Fig. |], together with the forward delay time 
from Eq. ( [2.19| ). The observed negative delay times for scattered waves reflect the negative 
derivatives of 5*-wave and other phase-shifts with respect to energy, weighted in the average 
time with contributions to the cross section. The falling phase-shifts make the interaction, 
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with regard to scattering, effectively repulsive above Ei a \, ~ 2.5 MeV. Notably, the NN 
interaction is often only considered repulsive when phase shifts are predominantly negative, 
although the energy derivatives of the phase shifts also necessitate a consideration. Both 
the T = and T = 1 phase shifts fall with energy, as the Levinson's theorem requires that 
the deuteron formation in T = channel and enhancement of the density of states in the 
d* region in T = 1 channel be compensated by a depletion in states at higher energies. In 
a semiclassical consideration, the delay time for scattering should be limited from below 
by —2d/v, where d is interaction range and v is relative velocity. While the negative delay 
times in Fig. f| decrease in magnitude with an increasing laboratory energy, their decrease 
is faster than implied by the above limit. At very low energies, the times for the scattered 
waves become positive. At Ei ab ~100 keV they begin to be governed by the singlet scattering 
length, At s ~ —a s /v. 

For comparison of angular dependence, the delay time for a hard-sphere repulsive- 
scattering (discussed more in the next section) as a function of the scattering angle is 



and averaged over forward and backward directions, given the constant scattering cross 
section, 



Correspondingly, for repulsive scattering, more negative delay times might be expected at 
6 = 90° than at 9 = 0. Indeed, that is observed in NN scattering at the higher laboratory 
energies. In fact, in the periphery the interaction might be expected attractive and the delay 
times might be expected to turn to zero or even positive in the forward directions, and not 
just decline in magnitude as for hard-sphere scattering, and this can be seen for the 360 MeV 
scattering in Fig. [j| In passing, let us note that to properly isolate forward and backward 
directions, constructing an amplitude prior to the antisymmetrization, one should continue 
the phase shifts over the missing partial waves. If one were to identify a hard-sphere radius 



At s (9) = sin -, 



(6.3) 



v 




(6.4) 
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d for simulations from the delay at 90°, quite low values would have been obtained compared 



range 10 - 30 MeV to d = 0.15 - 0.30 fm for E lab within the range 100 - 400 MeV. At low 
energies, the times become more negative at 0° than at 90°, and this is due to 3 Si- 3 Di 
interference. 

The forward delay times are positive in a wide energy range 3~^ a 6~150 MeV, see Fig. |j. 
In the lower portion of the range, this is due to the fact that, for large S'-wave phase-shifts, 
the real forward amplitude, multiplied by a relative momentum, increases, although the 
.S'-wave phase-shifts decrease. In the higher portion of the above energy range, when S'- 
wave contributions are low, the forward delay times are positive due to the positive energy- 
derivatives of phase shifts for some high partial waves which, unlike in the delays for scattered 
waves, are not weighted by partial cross sections in the forward direction. 

Of some interest is the issue of elastic NN interactions at very high energies when 
amplitudes are primarily diffractive |[L7| . Schematically, a purely diffractive amplitude may 
be represented, given Si = for £ < £ c and Si = 1 for £ > £ c , as 



As a diffractive amplitude is purely imaginary, its phase does not depend on energy and 
the time delays for the elastically scattered waves identically vanish, cf. Q2.14j ). Likewise, 
the times for the forward waves ( P-19| ) vanish. This is consistent with the concept of particles 
moving freely around the interaction region. 

We conclude this section with a presentation of the pressure in a low density nuclear 
matter at moderate temperatures, as a function of the density and temperature, such as 
should be, generally, reproduced in simulations. Within the second order in density the 
contributions to pressure (beyond the free-gas term P = n BT) come from the nucleon 
antisymmetrization, formation of deuterons, and nucleon-nucleon scattering, 



to what was used in simulations |fTE| , declining from d = 0.60 — 0.80 fm for within the 



f(0) = ^-J2(2£ + l)P e (cos9). 



(6.5) 




(6.6) 
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where the three respective contributions to the virial coefficient 02 (T) are given by 

2 5/2 a 2 (T) = 1-3 e^ r - £ (2T + l)(2j + 1) / ^ 2 ^ e"^, (6.7) 

TjtS 

and 5^ is the deuteron binding energy. The iViV antisymmetrization increases the pressure. 
Deuteron formation lowers the pressure, while scattering, at most temperatures, increases 
the pressure. The overall effect of interactions, following the Levinson's theorem, declines 
with the increase of temperature. The second virial coefficient ( |6.7| ) is shown as a function 
of temperature in Fig. [|. 



VII. DELAYS IN SIMULATIONS 

A. Effect of Scattering Prescriptions 

Without explicitly delaying or advancing particles as they pass near each other within 
a simulation, the time spent by particles in the vicinity of one another can be affected by 
the scattering prescription. The prescription can make At s positive or negative. Effects of 
scattering prescriptions on macroscopic features of reaction dynamics at Bevalac energies 
have been investigated in Refs. HT5HTB| , see also ||(J . We analyze two examples, hard sphere 



scattering and scattering at the point of closest approach. 

Our first example is that of hard-sphere two-particle scattering. Given a volume in 
relative coordinates of radius R, and a hard-sphere potential that rises as particles are 
distance d away, the expected reduction in the density of two-particle states within the 
relative volume, due to scattering, becomes Ap/p = At co11 /t = —d 3 /R 3 . Particles that 
make contact along a line at an angle a to the direction of original relative motion, get 
deflected by an angle 9 = n — 2a. When reaching the boundary of the relative volume 
defined with the radius R, the particles traverse a relative distance that is altered by the 
scattering. The alteration divided by the relative velocity, gives the change in the time spent 
in the vicinity of the other particle, Ar s (b), that depends on impact parameter and can be 
worked out from geometry, 
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The average change in time due to scattering is then given an integral over the product of 
the probability density that a specific scattering occurred, times the change in time in that 
scattering, 

a coil _ J dbbAr s (b) _ d 2 Ad 

Jdbb R 2 3v' 1 ' ' 

In the last expression, the factors d 2 /R 2 and (4/3) d/v represent, respectively, the probability 
that a collision occurs and the average time lost then in a collision. The average time spent 
within the volume of radius R follows from dividing the volume by cross-sectional area and 
velocity r = (4vr /3)R 3 /tcR 2 v = {A/3)R/v, and we find that At co11 /t = -d 3 /R 3 , as we 
expected. The pressure corresponding to (|673| ), ( |7.1| ) and (|7.2| ), in absence of statistical 
effects, from fl3.17| ) and ( |3.18| ), is 



P = nT + n 2 T — . (7.3) 
3 

Positive delay times are obtained when scattering the particles as if off a thin spherical shell 
of size d open in the direction of motion (case of a concave rather than convex mirror). Then 
expressions ( pM|) and (|7.2| ) remain valid but with changed signs. 



Within most common prescriptions for scattering in simulations, it is assumed that 



particles come abreast of each other while at a distance b < d = y a/ir. When deflected in 
a direction making an angle 9 relative to the original direction, within a plane at an azimuthal 
angle <fi with respect to the original plane containing particle trajectories, the particles reach 
a distance R^> d after a time longer by 

Ar s (6, St) = - - sin 9 cos , (7.4) 

v 

than in the absence of scattering. If the scattering is repulsive (0 = 0) and isotropic, 
i.e. cos ^-distribution is flat, then the average delay time for a given impact parameter 
becomes equal to 
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Ar.(&) = -|£. (7.5) 
By averaging over all impact parameters one further gets 

Notably, the time lost in any one collision is reduced here by only a factor 7r/4 compared to 
hard-sphere scattering. By choosing = n in collisions, one can produce the positive delay 
times. 

Generally, given required delay times such that |Ar s | < (2/3) d (sin 9) = (Ar s ) max , these 
times may be generated making a fraction v of all scatterings repulsive and a fraction \ — v 
attractive. This fraction is given by 

„ = i(l-*IL_ V (7.7) 



2 V (Ar s ) r 



B. Numerical Determination of Equations of State and Transport Coefficients 

While we have limited ourselves to the discussion of effects of two-particle interactions 
on thermodynamic properties, many-body calculations can provide information accounting 
for interactions of a few particles at a time |T0|JT2]] . The prescriptions for interactions in sim- 
ulations may affect microscopic thermodynamic quantities within a higher order than the 
second in density and they can also affect transport coefficients. Generally, given the pre- 
scriptions, thermodynamic quantities within simulations and coefficients can be determined 
numerically and confronted with those from fundamental calculations. 

All thermodynamic quantities can be derived once one knows the pressure as a function 
of the chemical potential and temperature, P(fi,T). To determine the pressure, the system 
may be enclosed in a box of macroscopic volume V, in contact with a free non-interacting gas 
with of chemical potential \x and temperature T, possibly only within some external potential 
lower than /i. The contact with the gas can be made through the walls in one direction, 
and in two remaining directions periodic conditions may be used. If clusters are produced 
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within the simulation, then the interfaces to the free gas can be made impermeable to those. 
Nucleons not in a cluster, on the other hand, getting into the free zone, would never return. 
At the same time, nucleons from the free zone with equilibrium phase-space distribution for 
given /i and T with the inclusion of statistics, would pour in into the box. The pressure 
within the box could be then computed using the following virial-type expression with terms 
for different possible ways of accounting for interactions in a simulation 

p = ~y \l ( / d tJ2 Pi • v i - S Pv ■ y ij At} 7 + J dt^ F lj ■ r y 

l\ j i<j i<j 

+ Apij ■ + / dt J dv (pu - J* {t) dp' uj I . 

The above equation is limited to the mean field being momentum independent. 



(7.8) 



With regard to (|7.8j ), pressure is, generally, recognized as the density of momentum flux 
in equilibrium, in any one direction. In (|7.8[) the pressure is evaluated by taking a trace 
of the momentum flux tensor and dividing it by 3 for the three directions; r is the time 
over which the system is investigated. The first term on the r.h.s. of ( |7.8j ), with a sum over 
particles in a box, accounts for the transport of momentum when particles move. The second 
term, with a sum over particle encounters, accounts for the situations when particles pass 
in the vicinity of one another and their relative motion is delayed in a simulation by At/. 
The vector pjj = (pj — Pj)/2 is relative momentum and Vjj = Vj — v.,- is relative velocity. 
The third term on the r.h.s. of ( f7.8|) accounts for transport of momentum in two-particle 
interactions treated explicitly. The force = is that due to particle j on particle i. 
Due to the interaction, the relative momentum changes over a distance r^. The fourth 
term accounts for instantaneous changes of relative momentum in collisions, and the final 
term accounts for the effects of interactions treated in the mean-field approximation. It is 
apparent in ( |7.8| ) that positive forward delay times, for a given particle number, reduce 
the pressure, and negative enhance. Further, the attractive scattering style, Apij ■ < 0, 
reduces the pressure, and repulsive, Ap^- • > 0, enhances. 

Transport coefficients such as shear viscosity or heat conductivity may be determined 
within a simulation by imposing different conditions within the free gas beyond the two 
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walls of the box in contact with that gas. Provided that the walls are perpendicular to the 
x-axis and at x — ±Ax, the viscosity coefficient may be e.g. determined by giving to the 
free gas the velocities in the ^-direction equal to ±Av, respectively, in the two separated 
regions. The coefficient then follows from an off-diagonal term of the momentum flux tensor 
in the box 

i<j i<j / 

To gain an insight into ( |7.9| ), one may consider a simple assessment of the viscosity in 
a medium when ignoring the effects of the finite range and duration of interactions, i.e. in- 
vestigating, in particular, only the effects associated with the first term in ( |7.9|) . Provided 
that particles propagate freely between collisions for an average time Tp ~ l/na(v), a parti- 
cle at a position Xi of velocity vf would have, on the average, a momentum in the y-direction 
such as characteristic for a position this particle had a time Tp earlier. With this, Eq. ( |7.9|) 
gives 

V ~ ~ VaI/Ax ? m {vV){Xl ~ V * Tf) V ' 

= ~ VaI/Ax ^ m ^~ < r ^ ( Av / Ax ) v ? 

// rx2\ mn(v 2 )T F m(v) ._ 

w mn ({v x Y) t f = N ' - » — — . (7.10) 

3 3cr 

The hard-sphere scattering is known to enhance viscosity to within the lowest order in 
density, but the enhancement factor for viscosity is smaller than that for pressure. For a given 
shear then the medium with hard-sphere scattering, to a lowest order, behaves as less viscous 
compared to one where effects of interaction range and duration may be ignored. For the 
numerical determination of a heat conduction coefficient, given assumptions of a simulation, 
the temperature should be set different in the two regions with free gas adjacent to the box 
where interactions take place. 
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VIII. CONCLUSIONS 



The principal goal of simulating heavy-ion collisions is to infer the equation of state of 
nuclear matter. The effective equation of state for a simulation depends on several aspects: 
treatment of the mean field in simulation, time delays for interactions and scattering pre- 
scriptions, inclusion of various resonances. In this paper we have carried out a detailed 
investigation of the time delays and of scattering prescriptions, and have shown that they 
should be incorporated in a coordination with the mean field. If that is not followed, the ef- 
fective equation of state may be inconsistent with two-body scattering which constrains 
thermodynamic quantities to within the second order in the virial expansion. 

In Sec. [n] expressions were derived for the average delay of an outgoing scattered wave as 



a function of the scattering angle. In Sec. pj, such delays in a dynamics were shown to be 
consistent with a two-body density of states if, besides, a forward delay or mean field were 
included and calculated in terms of the forward scattering amplitude. Other ergodically 
consistent prescriptions were presented where all the effective time delays were incorporated 
either entirely into the scattering or into the mean field. Alterations to these considerations 
for a Fermi-degenerate system were shown to be nontrivial in Sec. |V[ Sections |V] and |V] 
illustrated the time-delay considerations with the examples of resonance and NN scattering. 
In Sec. |V11| it was shown that repulsive and attractive scattering schemes can be interpreted 
in terms of time delays and equations of state associated with these schemes may be un- 
derstood quantitatively at a two-body level. The implications of scattering prescriptions for 
transport coefficients were also discussed and a practical method of determining the equation 
of state and coefficients for a simulation was presented. 

We conclude by giving some perspective to the considerations discussed here. Most of 
the flow in high-energy ~1 GeV/nucleon heavy-ion collisions stems from the one-body free 
gas pressure. The interactions of particles affect the pressure and can increase the observed 
flow of a collision by ~ 50% compared to a free gas. As experiments are now able to 
measure sidewards flow and squeeze-out flow differences to better than 20%, a detailed and 
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careful understanding of the simulations becomes crucial. Demonstrating that a simulation 
reproduces experimental results to within (10-20)% has little meaning unless thermodynamic 
properties of the simulation are understood to within (10-20)% as well. Simulations are 
becoming more sophisticated. Important aspects, such as e.g. the change of the dispersion 
relation for pions, are being incorporated through energy-dependent mean fields. At very 
relativistic energy, simulations must incorporate a large number of resonances which can be 
very broad compared to the temperature. It is hoped that the prescriptions and constraints 
presented here will contribute to both the development of improved codes and to a better 
understanding of existing approaches. 
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APPENDIX: 

In this Appendix we investigate changes in two-particle scattering-rate induced by 
the presence of other particles, within the lowest order in density, in the context of er- 
godic theorem. To reach directly essential results we adopt simplifying assumptions. First, 
we concentrate on a subsystem of one of the particles engaged in scattering, denoted as 1, 
and the spectator particle denoted as 2. We assume that these two particles interact only 
perturbatively with all other particles within the system, but not necessarily with one an- 
other. If the system has many particles confined to a finite volume and sampling is done 
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over a long time r, then we may write an ergodic theorem for the subsystem of particles 1 
and 2 alone: 

1 dr ( r dEl dn^ EI/T Y l J^ e -E/T_ (A1) 



r dp dE \J dE' J dp dE 

In the above dr/dpdE denotes time spent in a given set of internal states, per unit volume 
of energy-momentum space, dn/dpdE denotes the number of states, and e~ E ^ T is relative 
probability that any single internal state is occupied at any instant. An equivalent formula- 
tion of the ergodic theorem (|A1|) is that for any two states (p, E) and (p', E')\ 

( dr \ 1 dn WIT ( dr \ 1 dn p , /T , A . 

«-^= — — - —— e - E f T . (A2) 



\dpdEJ dpdE \dp' dE' J dp' dE' 

For simplification, we next assume that 2 is much heavier than 1, whereupon we can make no 
distinction between relative states of 1 and 2 and single-particle states of 1. Equation ( |A2|) 
then may be also written as 



( dT ) 1 dRl e'Ei/T = ( dT ) 1 dn l e -E>/T (A3) 

\dp 1 dE l j dp 1 dE 1 {dp^dE'J dp' l dE' 1 v ; 

We are now set to address the scattering of a particle 3, representing remainder of 
the system, with particle 1, in the presence of 2. Transition or scattering rate within the 
system is, generally, given by a transition matrix element squared multiplied by the density 
of final states in energy, times 2tc. When sampling is carried out over a long time, then 
the number of transitions within the system from one set of states to some other should 
be the same as from the other set to the first. On multiplying both sides of Eq. (|A3|) by 
common factors and after manipulations, we can demonstrate an equivalence of the ergodic 
condition (|A3|) with the general condition of the equality of the number of transitions in 



the different directions. Specifically, we multiply both sides of ( |A3| ) by a product of average 
number densities of 3, dN 3 /dp 3 dE 3 = dn 3 /dp^ dE 3 x e~ E3 ^ T , in the vicinity of two different 
states (p3, E 3 ) and (pg, E' 3 ), such that pi + p 3 = + P3 and Ei + E 3 = E[ + E' 3 , and by 
a matrix element squared in the momentum representation, same for direct and inverse 
transitions, |A^i3_»i'3'| 2 = 1 2 = \Ai\ 2 . After manipulations we get 
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dr dN 3 (27r) 4 . |2 dD 3 



dp! dEi dp 3 dE 3 V 3 1 1 dp[ dE[ dp' 3 dE' 3 

dr dN 3 (2?r) 4 2 dni dn 3 
dpi dE[ dp' 3 dE' 3 V 3 1 1 dp x dE x dp 3 dE 3 ' 1 ' 

The volume V is assumed here to be large compared to the range of interactions. 
The l.h.s. of (|A4|) represents number of transitions that take place, per element of energy- 
momentum space dpi dE\ dp 3 dE 3 dp[ dE[, from states 13 to 1'3'. The r.h.s. represents num- 
ber of inverse transitions per energy- momentum element dp[ dE[ dp' 3 dE 3 dpi dE\. By virtue 
of energy-momentum conservation the energy-momentum elements involving any three out 
of the four states in the above are actually identical. The number of transitions in ( |A4[ ), 
e.g. on the l.h.s., is represented in terms of the time spent by particle 1 in the region (p 1; E\) 
times the probability of finding particle 3 in the region (p 3 ,E 3 ), times the transition rate 
into (p'^Ei) and (p' 3 ,E' 3 ). 

Equation (|A4|) is fully equivalent to ( |A3| ) or (|A1|) . Given that interaction of particle 3 
with any other particles is perturbative, the perturbative scattering of 1 and 3 is modified 



in the presence of 2, compared to free space, by a changed final density of states. From (|A4|) 
it follows that, to comply with ergodicity, it is necessary to allow particles to participate 
in transitions at any time when in a given relative state. In the transition rate and, in con- 
sequence, in the cross section, it is necessary to account for the changed final-state density. 

Within many-body theory, a final-state density for scattering is typically described in 
terms of single-particle spectral functions [[H],^] equal, up to a factor, to the imaginary 
part of single-particle Green's function ( |3.8| ), A = — 21m g. For example, for the case 
above, the number of transitions in which particles from (pi,Ei) and (ps,^) interact and 
populate (p[,E[) and (p 3 ,E 3 ), per unit time and per momentum volume dp 3 dp[, given 
well-defined energies of 3, would be represented as 

j^y 6 f(P3)\M\ 2 A(p' 1 ,E[). (A5) 

On comparing the previous expression for transitions with the one in terms of A, one can 
conclude that, in the discussed case, 
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V N dni dn 

T^ A ^ El) = ^LEr^m- (A6) 

This implies an integral relation between the density of relative states p and A, 

dn V f 

wr p ^ ) = W)il d ^ A ^ E ^ (A7 » 

We now proceed to examine an explicit form of A in the T-matrix approximation in which 
interactions within a two-particle system, such as that of 1 and 2 in the above, are fully 
accounted for. With an explicit form of A, we shall further examine the validity of ( |A7| ), 
establishing a correspondence between different terms in A and in p. We shall then discuss 
in-medium scattering rates. 

From (13.81) we find 



A( Pl , E x ) = -2 lmg( Pl , E 1 ) = — -l^^l— ___ . (A8) 

(Ei - e(pi) - u(pi, Ei)) 2 + (7(pi, Ei)) 2 /4 

Here the real and imaginary parts of self-energy separately depend on energy and momentum. 
In the T-matrix approximation, this dependence corresponds to a separate dependence of the 
T-matrix on energy and momentum as in the Lippman-Schwinger equation. On expanding A 
for low scattering rate 7, in the discussed case rate for scattering of 1 off 2 given large V, 
we obtain 



A(p 1 ,E 1 ) 



-1 

du(p 1 ,E 1 



dE 1 



2tt5(E 1 -£(p 1 )) 



Ei=£(pi), 

V 



- 7 (pi,Ei) — (A9a) 
Ei — t{Pl) 

9u(pi,Si) 



2 7 r5(E 1 -£(p l )) 



2tt5(E 1 -S( Pi )) 

Bi=£(pi) 



- f(p " £i) ir^)- (A9b) 

The energy £(pi) in the above is a solution of the equation E\ — e(pi) — u(pi,Ei) = 0. 
The factor multiplying the 2ir 5 in ( |A9a| ) is termed a wave-function renormalization factor, 
and 
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V d V d 1 / 1 1 \ d 1 , . . , 
— = = — hm - H = hm — — : . (A10) 



x dx x dx e->o 2 \x + ie x — ie) e->o de 2i \x + ie x — ze 



see also |4T| , f42"fl . Equation ( |A9b| ) follows from ( |A9a| ) on recognizing that, according to a dis- 
persion relation, the energy derivative of the mean field is proportional to the off-shell 
scattering rate. Thus, it is expected small when the rate is small. If V is large enough, then 
the mean field itself is small, and we can further expand A in terms of u, 

A(p u E x ) « 2tt 5(E 1 - e(pO) - u(p ls e( Pl )) 2vr 5'(E 1 - e( Pl )) 
du(p 1 ,E 1 



+ 



2n5(E 1 - e(p 1 ))- 1 (p 1 ,E 1 ) — (Alia) 

E 1= e( Pl ) E l - 



dE x 

= 2tt {8{E x - e(p x )) - u( Pl , E x ) 5'iEt - e{p x ))) 

-7(pi,^i) w v (AUb) 

^i - e(pi) 

Within the T-matrix approximation, the mean field and scattering rate due to 2 are 

u( Pl ,J50 = / ^/(p 2 )He< P |r(J5)|p> = He< Pl |r(J5 l )|p 1 >/^ (A12) 

and 

7 (pi,^i) = -2/ ^/( P2 )Im(p|T(£)|p) = -2Im(p 1 |T(^ 1 )|p 1 )/F, (A13) 

with center expressions following under assumption of well-defined energies of 2, and r.h.s. ex- 
pressions representing results for the specific case under discussion. 

On inserting A in the form ( |Alla| ) into the r.h.s. of ( |A7|) , we find that the leading term 
in (|Alla|) produces po from Eq. ( |3.1|) . Second term in (|Alla| ), with a derivative of the 



5-function, gives a complete contribution to Ap associated with the forward delay time, 

^U^(p\ReT(E)\p^ =i£(2£ + l)cos2^, (A14) 



compare Eq. (|3.11| ). The derivative from the 5-function acting on po gives a contribution 
associated with ATf as , cf. Eq. Q3.4| ), and the derivative acting on u gives a contribution asso- 
ciated with Ar 2 das , cf. Eq. (gT5j) . Finally, the remaining wave-function renormalization and 
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off-shell terms in ( |Alla|) yield jointly, on insertion into the r.h.s. of (|A7|) and after lengthy 
manipulations, a scattering contribution to p of the form p av Ar s /V, where At s is given 
by ( |3.12| ), compare ( |3.2| ) and (|3.7|), indeed confirming the equality in (|A7|). The manipula- 
tions involve, in particular, expressing the <5-function and principal value in ( |A1 1| ) in terms 
of the imaginary and real parts of the free two-particle Green's function, and an extensive 
use of the relations between the imaginary and real parts of the free two-particle Green's 
function and the T-matrix, following from the Lippman-Schwinger equation, Eq. (|4.8| ) with 
AT = 1 in Go- Presence of the momenta pi in the last term in (|Alla| ) and in ( |A13| ), far 



from the shell defined by e(pi) = E\, indicate the effects of near zone in the interaction of 1 
and 2. 

Generally, in a scattering rate a single-particle spectral-function A would be used for 
each of final-state particles. On the basis of the example above, one would conclude that 
use of such function leads to the population of states consistent with two-particle density. 
A more general relation between A and p, than ( |A7| ), is 

27 / dp JjLz) A ( Pl ' E + P2/2M ~ e °° 2) ) = Ap{E)} (A15) 
where integration is carried out over relative momentum at fixed total momentum P = 
Pi + P2- Use of the distribution /(P2) in (|A12j) and ( |A13| ) allows for various momentum 



values of 2. As / is not normalized to yield one particle within a given volume, there is 
a possibility for 2 being absent from V, giving a reduction in the relative weight of the 



correction to the density of states due to the interaction, cf. (|A11 ) 



In practical applications, the incorporation of the first correction term in A in (|Alla|) 



associated with the forward delay time, amounts to correcting single-particle energies in 
the scattering rates by the mean field. The incorporation of the scattering delays in the 
final states in scattering rates can be much more cumbersome given the form of terms 
in ( |Alla|) . One possible solution to get consistency with ergodicity is to multiply rates or 
cross sections by factors of the form 

1 + J dp 2 f(p 2 )avAr s , (A16) 
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for each of final-state particles. That would be analogous to putting all delays into the 
forward delay time in the particle propagation, discussed in Sec. [] fD|. A weakness of an 
approach where just final-state densities in two-particle scattering are modified on account 
of scattering with other particles, is the disregard of correlations that may persist through- 
out the interaction process. Effects of correlations become even apparent in a more detailed 
analysis of the discussed simple example with perturbative scattering. Thus a careful reader 
might notice that the matrix elements in flAJ ) for transitions should not be, generally, taken 
between plane waves but rather between eigenstates of an internal hamiltonian of 1 and 2. 
Separation of the wave-functions for these states into incident and scattered portions would, 
generally, yield two-particle and three-particle scattering terms in the transition rate. In- 
terference of the forward and scattered waves for 1 and 2 would yield shadowing corrections 
in the matrix element for the scattering of 1 and 3. The genuine three-particle scattering 
term would be characterized by the lack of momentum conservation within the subsystem 
of 1 and 3 alone. Overall, such effects are, though, beyond the scope of present paper. 
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FIGURES 

FIG. 1. Time delays in the case of a Breit-Wigner resonance as a function of energy from 
the resonance divided by the width. The solid line represents the time for the scattered wave 
Ar s . The short-dashed line represents the time At' s which is ergodically consistent when the delay 
for forward wave or mean field are neglected. Finally, the long-dashed line shows the time delay 
for the forward wave divided by the fraction of the incoming wave that is scattered, At' s — At s 
(cf. Eq. 0). 

FIG. 2. Time delays for a ir + p system as a function of cm. kinetic energy. The solid line 
represents the time delay for the scattered wave averaged over angles and spin directions, At s . 
The dashed line represents the forward time delay averaged over spin directions and divided by 
the fraction of the incoming wave that is, on the average, scattered, At' s — At s . 

FIG. 3. Spin-isospin averaged time delay for NN scattering as a function of the cosine of 
cm. scattering angle, at several values of laboratory kinetic energy indicated in the figure in MeV. 

FIG. 4. Time delays for NN system as a function of laboratory kinetic energy. The solid line 
represents the time for the scattered wave At s averaged over the scattering angle and the spin and 
isospin directions. The dashed line represents the time for the forward wave, averaged over spin 
and isospin directions and divided by the fraction of incident wave scattered, At' s — At s . 

FIG. 5. Second virial coefficient for nuclear matter, multiplied by 2 5 / 2 , from Eq. (|6.7| ), as 
a function of temperature. 
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